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Abstract. A deformed relativistic Hartree-Bogoliubov theory in continuum has been developed for 
the study of neutron halos in deformed nuclei and the halo phenomenon in deformed weakly bound 
nuclei is investigated. Magnesium and neon isotopes are studied and some results are presented for 
the deformed neutron-rich and weakly bound nuclei “+Mg and *°Ne. The core of the former nucleus 
is prolate, but the halo has a slightly oblate shape. This indicates a decoupling of the halo orbitals 
from the deformation of the core. The generic conditions for the existence of halos in deformed 
nuclei and for the occurrence of this decoupling effect are discussed. 
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INTRODUCTION 


The halo phenomenon is one of the most interesting topics in modern nuclear physics. 
Since most open shell nuclei are deformed, the interplay between the deformation and 
the weak binding feature raises interesting questions. In order to give an adequate 
description of possible halo in a deformed nucleus, a model should be used which 
includes in a self-consistent way the continuum, large spatial distributions, deformation 
effects, and couplings among all these features. For this purpose, a deformed relativistic 
Hartree-Bogoliubov theory in continuum has been developed and the halo phenomenon 
in deformed weakly bound nuclei is investigated. 


THEORETICAL FRAMEWORK 


Continuum contribution and relativistic Hartree-Bogoliubov theory 


By solving the non-relativistic Hartree-Fock-Bogoliubov (HFB) [1, 2, 3] or the rela- 
tivistic Hartree Bogoliubov (RHB) [4, 5, 6, 7, 8] equations in coordinate (r) space, the 
mean field effects of the coupling to the continuum can be fully taken into account. With 


the relativistic continuum Hartree-Bogoliubov (RCHB) theory [7, 8], properties of the 
halo nucleus !!Li has been reproduced quite well [4] and the prediction of giant halos 
in light and medium-heavy nuclei was made [9, 10, 11]. The RCHB theory has been 
generalized to treat the odd particle system [12] and combined with the Glauber model, 
the charge-changing cross sections for C, N, O and F isotopes on a carbon target have 
been reproduced well [13]. 

The Dirac Hartree Bogoliubov (RHB) equation for the nucleons reads [14], 
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where E, is the quasiparticle energy, À is the chemical potential, and Ap is the Dirac 
Hamiltonian, 
hp(r,r’') =a@-p4+V(r)+B(M+S(r)). (2) 


with scalar and vector potentials 
S(r) a 8go0(r), (3) 
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The equations of motion for the mesons and the photon 
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have as sources the various densities 
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where, according to the no-sea approximation, the sum over k > 0 runs over the quasi- 
particle states corresponding to single particle energies in and above the Fermi sea. 
In the particle-particle (pp) channel, we use a density dependent zero range force, 


VPP (ri, ro) = Vox (I~ P°)8(r1 — r2) (1-2) : (7) 


5(1 — P?) projects onto spin S = 0 component in the pairing field. The pairing potential 
then reads, 
A(r) = Vo(1 — p(r)/Psat) K(1), (8) 


and we need only the local part of the pairing tensor 
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Large spatial density distribution and the Woods-Saxon basis 


In order to (1) consider properly the asymptotic behavior of nuclear densities at 
large r and (2) make the numerical procedure less complicated, the Woods-Saxon basis 
has been proposed in Ref. [15] as a reconciler between the harmonic oscillator basis 
and the integration in coordinate space. Woods-Saxon wave functions have a much 
more realistic asymptotic behavior at large r than do the harmonic oscillator wave 
functions. On one hand, one can still use a large box boundary condition to discretize the 
continuum and easily find numerical solutions for a spherical Woods-Saxon potential 
in r space. One the other hand, these Woods-Saxon wave functions can thus be used 
as a complete basis for spherical or deformed systems, and one finally comes back to 
the familiar matrix diagonalization problem. It has been shown that the Woods-Saxon 
basis can satisfactorily reproduce the large neutron density distribution in weakly bound 
nuclei obtained in r space [15]. Recently, for spherical systems, both non relativistic 
and relativistic Hartree-Fock-Bogoliubov theories with forces of finite range have been 
developed in a Woods-Saxon basis [16, 17]. 


Deformed relativistic Hartree Bogoliubov theory in continuum 


The Woods-Saxon basis can be extended to more complicated situations for exotic 
nuclei where both deformation and pairing have to be taken into account. Over the 
past years, lots of efforts have been made to develop a deformed relativistic Hartree 
theory [18] and a deformed relativistic Hartree Bogoliubov theory in continuum (the 
DefRHBC theory) [19, 20, 21, 22, 23]. 

For axially deformed nuclei with the spatial reflection symmetry, we expand the 
potentials S(r) and V(r) in Eq. (2) and various densities in terms of the Legendre 
polynomials [24], 


=È flr r)P,(cos@), A = 0,2,4,- (10) 
with an explicit definition of f} (r). 


The quasiparticle wave functions Ug and V; in Eq. (1) are expanded in the Woods- 
Saxon basis [15]: 


Uc(rsp) = Yup’? Paxm(rsp); (11) 
nk 


Vi(rsp) = Loi a Paxm( TSP). (12) 
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Onxm(rsp) is the time reversal state of @:xm(rsp). Because of the axial symmetry the 
z-component m of the angular momentum j is a conserved quantum number and the 
RHB Hamiltonian can be decomposed into blocks characterized by m and parity 7. For 
each m”-block, solving the RHB equation (1) is equivalent to the diagonalization of the 
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and 
oe (HDi) = ((nkm|hp|n'x’,m)), (15) 
B = (Ahu) = (amA, m). (16) 


Further details are given in the appendixes of Ref. [23]. 

In order to describe the exotic nuclear structure in unstable odd-A or odd-odd nuclei, 
the DefRHBC theory has been extended to incorporate the blocking effect due to the 
odd nucleon(s) [25] The deformed relativistic Hartree-Bogoliubov theory in continuum 
with the density-dependent meson-nucleon couplings is developed recently [26]. 


RESULTS AND DISCUSSIONS 


We next present some results from the DefRHBC theory by taking magnesium and neon 
isotopes as examples and discuss some results for the deformed neutron-rich and weakly 
bound nuclei “Mg and *°Ne [21, 22, 23]. 

Magnesium isotopes have been studied extensively in Refs. [21, 22, 23] with the 
deformed relativistic Hartree-Bogoliubov theory in continuum and the parameter sets 
NL3 [27] and PK1 [28]. For the pp interaction (7), the following parameters are used: 
Psat = 0.152 fm~> and Vo = 380 MeV-fm?, and a cut-off energy EJP ` = 60 MeV is 
applied in the quasi-particle space. These parameters were fixed by reproducing the 
proton pairing energy of the spherical nucleus 7°Mg obtained from a spherical relativistic 
Hartree-Bogoliubov calculation with the Gogny force D1S. A spherical box of the size 
Rmax = 20 fm and the mesh size Ar = 0.1 fm are used for generating the spherical Dirac 
Woods-Saxon basis [15] which consists of states with j < oh. An energy cutoff Ex), = 
100 MeV is applied to truncate the positive energy states in the Woods-Saxon basis and 
the number of negative energy states in the Dirac sea is taken to be the same as that of 
positive energy states in each (£, j)-block. 

In our deformed RHB calculations with the parameter set NL3, 4°Mg is the last 
nucleus of which the neutron Fermi surface is negative and the two neutron separation 
energy is positive [21]. In the calculations based on the parameter set PK1, *7Mg is the 
last bound nucleus in Mg isotopes [23]. 

It was found in Ref. [21] that the ground state of “Mig is well deformed with 
quadrupole deformation 2 = 0.32 and a very small two neutron separation energy 
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FIGURE 1. (Color online) Density distributions of “Mg with the z-axis as symmetry axis: (a) the 
proton density (for x < 0) and the neutron density (for x > 0), (b) the density of the neutron core, and (c) 
the density of the neutron halo. In each plot, a dotted circle is drawn for guiding the eye. This figure is 
originally published in Ref. [21]. 


San = 0.44 MeV. In the tail part, the neutron density extends more along the direction 
perpendicular to the symmetry axis. The density distribution is decomposed into contri- 
butions of the oblate “halo” and of the prolate “core” as shown in Fig. 1. The density 
distribution of this weakly bound nucleus has a very long tail in the direction perpendic- 
ular to the symmetry axis which indicates the prolate nucleus “Mg has an oblate halo 
and there is a decoupling between the deformations of the core and the halo. 

As discussed in Refs. [21, 22, 23], the shape of the halo originates from the intrinsic 
structure of the weakly bound or continuum orbitals. It turns out that in “4Mg, the 
essential level of the halo has a large contribution from the prolate A = 1 (p wave) 
component. In Ref. [21] an example was also discussed in which the halo and the core 
have similar shapes. In Fig. 2a the density distributions of all protons and all neutrons in 
the prolate deformed nucleus *°Ne are shown (b2 = 0.52). It can be seen that the neutron 
density not only extends much farther in space but it also shows a halo structure. The 
neutron density is decomposed into the contribution of the core in Fig. 2b and that of the 
halo in Fig. 2c. In contrary to the nucleus “Mig, we observe now a prolate halo, because 
the essential level of the halo has a large contribution from the prolate A = 0 (p wave) 
component. 

The halo feature is connected with relatively large cross sections and narrow lon- 
gitudinal momentum distributions in knockout reactions. The decoupling between the 
deformations of the core and the halo may manifest itself by some new experimen- 
tal observables, e.g., the double-hump shape of longitudinal momentum distribution in 
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FIGURE 2. (Color online) Density distributions of *°Ne. Details are given in Fig. 1. This figure is 
originally published in Ref. [21]. 


single-particle removal reactions and new dipole modes, etc. In particular, a combina- 
tion of the experimental method proposed in Ref. [29] and the theoretical approach de- 
veloped in Ref. [30] would be useful in the study of longitudinal momentum distribution 
in single-particle removal reactions with deformed halo nuclei as projectiles. The shape 
decoupling effects may also has some influence on the sub-barrier capture process in 
heavy ion collisions [31]. 

For odd particle system, the formation and the size of a halo depend on the interplay 
among the odd-even effects, continuum and pairing effects, deformation effects, etc. 
Some progress on this topic has been made recently [32]. 


SUMMARY 


We present recent progresses of the development of a deformed relativistic Hartree- 
Bogoliubov theory in continuum (DefRHBC) and the study of neutron halo in deformed 
nuclei. In the very neutron-rich deformed nucleus “Mg, pronounced deformed neutron 
halo was found. The halo is formed by several orbitals close to the threshold. These 
orbitals have large components of low ¢-values and feel therefore only a small centrifugal 
barrier. Although “+Mg and its cores is prolately deformed, the deformation of the halo 
is slightly oblate. This implies a decoupling between the shapes of the core and the 
halo. The mechanism is investigated and it was concluded that the existence and the 
deformation of a possible neutron halo depends essentially on the quantum numbers of 
the main components of the single particle orbits in the vicinity of the Fermi surface. 
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